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Abstract 

We study primordial tensor power-spectra generated during inflation in bimetric gravity. More 
precisely, we examine a homogeneous expanding spacetime in a minimal bimetric model with 
an inflaton and calculate tensor pertnrbations on the homogeneous background under slow-roll 
approximation. In terms of the mass eigenstates, only the power-spectrum of the massless state 
remains constant and both the power-spectrum of the massive state and the cross power-spectrnm 
rapidly decay dnring inflation. The amplitude of the physical power-spectrum is suppressed due 
to the flavor mixing. All power-spectra in the flavor eigenstates coincide with each other up to the 
first order of the slow-roll parameter. 
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I. INTRODUCTION 


Although no one doubts the existence of gravitons, no one has really observed them. 
Indeed, we do not know if gravitons have mass or how many species there are. What we 
can conhdently say is at least one of them must be sufficiently light in order to realize the 
Newtonian potential. Thus, we have room to suppose two or more graviton species exist. It 
is well known that two interacting massless graviton can not exist. However, it was recen^ 
found that a massless graviton and a massive graviton can exist at the same time liJJ. 
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The form of interaction terms is highly constrained in order to avoid ghosts. Such theory 
necessarily includes another metric in addition to the physical metric and is called bimetric 
gravity. In general, the interaction terms include hve theoretical parameters. 

We are interested in if the bimetric theory is theoretically consistent and reconciles with 
known experiments. One of important check points is if the predictions coming from inflation 
in bimetric gravity are consistent with cosmological observations which are now getting more 
precise. We examined the homogeneous expanding solutions and their stability in the slow- 
roll limit in our previous paper jol. We studied minimal bimetric models and obtained 
the unique stable branch of the solutions. Now, we construct homogeneous inflationary 
solutions under slow-roll approximation which correspond to the stable branch. Then we 
consider tensor perturbations on the homogeneous solutions and we calculate primordial 


;he slow-roll parameter. 

, Isj, treating perturbations on 
12| and on the detectability of gravitational waves 13| in bimetric 


tensor spectra generated during inflation to the first order o 
We note there are works dealing background solutions 


the FLRW background jo 
gravity. 

This paper is organized as follows. In section HIl we introduce an inflaton to a minimal 
bimetric model and we construct inflationary background solutions. We explain the prop¬ 
erties of the functions which are specific to bimetric gravity. We briefly mention slow-roll 
approximation and introduce a slow-roll parameter we use in the following calculation. In 
section m we derive the second order Lagrangian for tensor perturbations. There are two 
views of this system: the flavor eigenstates and the mass eigenstates, and we use the mass 
eigenstates for simplifying calculations. In section IIVI we calculate the primordial tensor 
spectra up to the first order of the slow-roll parameter by making use of the interaction 
picture. In the hnal section, we discuss the features of the tensor power-spectra in bimetric 
theory. 


II. INFLATIONARY BACKGROUND SOLUTIONS IN BIMETRIC GRAVITY 

In this section, we construct inflationary background solutions with a scalar held coupled 
to the physical metric in bimetric gravity. We write a minimal bimetric action with a 
canonical scalar held. We substitute a homogeneous isotropic metric ansatz into the action 
and derive the equations of motion by using the variational principle. We mention the 


3 










features of the solutions based on our analysis in our previous paper 




A. Action and Ansatz 


We consider a minimal bimetric action including a scalar field coupled to the physical 
metric and we substitute a homogeneous ansatz for metrics and the scalar held into the 
action. 


We use as the physical metrics 


fnu as the other metric and as the scalar held. 


fiu 

When we consider a bimetric theory 1^, the form of the interaction terms are restricted 
in order to avoid the Boulware-Deser ghost which include hve theoretical parameters 
{an\{n = 0,1, • • • ,4). For simplicity, we set 02 = 1 and other four theoretical parameters 
are equal to zero. Then we have 
M 2 


iVj ^ C P 1 




1/ fli 


( 1 ) 


where Mg and Mj are the Planck scales of the physical metric and the other metric respec¬ 
tively and we dehned the reduced Planck scale as 


I/Ml := 1/M| + 1/Mj . 


( 2 ) 


Indices written in Greek letters run over 0, • • • ,3. The last part in the action includes the 
interaction terms of the physical metric and the other metric where we dehned 


K ■= K - . 

We impose the homogeneous isotropic ansatz for these metrics as 

g^udx^dx’' = —N'^{t)dt^ + e^^^^'^'jijdx^dx^ , 


( 3 ) 


( 4 ) 


fgvdx^dx'^ = —M‘^{t)dt^ -f e^^^^^'jijdx^dx^ , 


( 5 ) 


where 7 *^ is three-dimensional hat metric and indices written in Roman letters run over 
spatial coordinates, i.e. z = 1, • • • , 3. N and M are lapse functions and e“ and are scale 
factors which depend only on time. We require the scalar held is also homogeneous 


ip = p{t) 


( 6 ) 
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When we substitute the ansatz into the Lagrangian, we hud 


L 


2„3a 


Mte 


3q! 


A,2 


N 


+ e 


3a 


2N 


- NV{ip) 


+ 


3/32 

'Id 


[iV(6 - 9e + 3e2) + M(-3 + 3e)] , 


( 7 ) 


where dots denote the time derivative and e is the ratio of the scale factors of the physical 
metric and the other metric, i.e. 


e:= . (8) 

We can see that a, 13 and (p are dynamical variables and N and M are non-dynamical 
variables included in the Lagrangian linearly. 


B. Equations of motion 


We derive the equations of motion of the dynamical variables a, f3 and p and two con¬ 
straints from the variational principle. We have a relation between the lapse functions 
in order for the two constraints to hold during time evolution 15j. We also describe the 
behavior of the solutions of the equations. 

From the variations of the Lagrangian with respect to the dynamical variables a, {3 and 
<p, we obtain 

1 


lf“V = 

N\N 


2 2/3 \ / M 

m cos 2 


2M2 \NJ 




(9) 


m\m) 


sin^ a( - — e ) (1 


_ 

~ 'm)'^ ’ 


( 10 ) 


N\n) ^ dp 


( 11 ) 


respectively, where we dehned a new parameter tana := Mg/Mf. Since N and M are non- 
dynamical and included only linearly in the Lagrangian, we obtain two constraints from the 
variations of the action with respect to them. One of them is 


tt \ 2 0 0/ On 1 


nJ 


) = 77/2 cos2 a(—2 -|- 3e — + 


3M| 




( 12 ) 


5 















which comes from the variation of the action with respect to N. The variation of the action 
with respect to M yields the other one 


2 . 2 1 
= m sm a— 


M 


( 13 ) 


We set the time derivative of eq. flT^ is equal to zero so that the constraint is satisfied during 
time evolution. We combine it with eq.(l9j) and obtain the following equation 


=0 • 


We can obtain the same equation also by using eq. flTdl) and 


eg 


is equal to zero, the solutions are known to be pathological [16 
the second factor is equal to zero, 

dd 

M = -^eN = (eN 
da 


where we defined 


C:= 


da 


( 14 ) 

Toll . When the first factor 
20|. Therefore, we assume 

( 16 ) 


( 16 ) 


Since iV is a gauge variable, we can set an arbitrary value for N. Substituting this relation 
into eq. (lT3|) leads to 

1-e 


d \ 2 

n) 


1 = sin^ a- 


( 17 ) 


This is another expression for the Hubble expansion of the physical metric. We can see that 
the Hubble expansion is a function of e. By equating eq. llT^ and eq. llTTl) . we obtain an 
equation 


sin^ a 


1-e 


= cos^ a(—+ 3e — 2) + 


1 


3M| 


1 143' 2 




( 18 ) 


This equation says that we find the values of e if we determine the energy density on the 
physical spacetime. We have examined the properties of the roots of this equation in the 
slow-roll limit in our previous paper jo]. In the slow-roll limit, we neglect (p contribution and 
this equation reduces to algebraic equation with constant coefficients, therefore, the roots 
are constant. This equation basically has three roots since it is a cubic equation. One of the 
roots is always negative, therefore it is not appropriate solution when we take into account 
the definition of e. Another is always larger than 1, therefore it is also inappropriate because 
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the Hubble expansion, which is written as eq. ffT7|) . becomes imaginary. The other one has 
the value between 0 and 1, therefore it is the only adoptable root as de Sitter spacetime. 
We have conhrmed that the root satishes the Higuchi bound 2l|, l22|, which is the stability 
condition of de Sitter spacetime with a massive graviton. 

Substituting eq. (fT^ into eq.([n]) and eq. lfllJ]) . we obtain the equation determining the value 

of C, 




C = id 

where we dehned the effective mass 

mlii(e) := rrf 


V|K»(£)-2(f)2] 


e cos^ a H— sin^ a (3 — 2e) . 


(19) 


( 20 ) 


In the slow-roll limit, C, is equal to 1 and therefore a and 13 are different only by a constant. 
Using the dehnition of we have a relation between ( and e as 


5C := C - 1 = - • 

ae 

We can see that vanishes in the slow-roll limit. 


( 21 ) 


C. Slow-roll approximation 


We introduce a slow-roll parameter and explain the assumptions used in the following 
analysis. We use the iV = 1 gauge in this subsection. 

We dehne a slow-roll parameter as 


s : = 


H 

IP ■ 


where H := a. From eq. ffT^ . we hnd 

<5C = 




2sH^ 


~ 2hf^) rnlg — 2 iJ 2 _ rn‘^{cos‘^ a)e{3 — 2e) 


= 2s(l -e) 


( 22 ) 


(23) 


The third expression is obtained by eliminating 0 with eq.(l9]) and the hnal expression is 
obtained by using eq. (fT7D and eq. fl20|) . From this equation, we can read off an expression of 
if as 


= 2s(l - - 2H^) = 2M^i-H){l -e)^j^-2)= 2M^{-H) (1 + 


m 


6^(3 -2e) 


tan^ a 


(24) 
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The third expression is obtained by using the definition of s and the last expression is 
obtained by using eq. ffT7|) and eq. fl20|) . We can see that the expression for ip is slightly 
changed from the conventional case . 

We assume the slow-roll parameter is small and constant. We neglect the higher order 
contribution O(s^) in the following. Under this approximation, we can see 5C, is constant 
i.e. C is constant but e has time dependence. We can easily obtain the explicit form of e by 
integrating the differential equation e = 2seH{l — e) obtained from eq. fl2T]) and eq. fl2^ : 


( 25 ) 


(1 - eo) + 


eo is the value of e in the slow-roll limit. We note that the functions with the subscript 0 


have the values in the slow-roll limit and all of them are constant in the following. 

III. SECOND ORDER LAGRANGIAN FOR TENSOR PERTURBATIONS 


We derive the second order Lagrangian for tensor perturbations on the homogeneous 


isotropic inflationary background. First, we obtain the Lagrangian in the flavor eigenstates, 
and then, we move to the mass eigenstates. 

A. Flavor eigenstates 

We give tensor perturbations to the original Lagrangian and derive the second order 
Lagrangian in the flavor eigenstates. We will see the variables in the flavor eigenstates are 
not decoupled from each other in the slow-roll limit. 

We consider perturbations such as 


Sgij = , Sfij = 


( 26 ) 


which satisfy the following transverse traceless conditions: 



( 27 ) 


We can decompose the tensor perturbations using polarization tensors as 


Qij = Qieij , Pij = Pielj , 
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(28) 



where the subscripts / correspond to the plus mode and the cross mode, i.e. / = {+, x}. 
When we substitute the tensor perturbations into the action, the second order Lagrangian 
is reduced to 


^ 4 

/ 


12 m 


1 




1 


--Ne^^ 
4 

4 


1 Pi 
L2M2 

1 / 99 ^ 2 

L2 ViV 


2e2“ 

2 1 e 


q! + 


■ 0/ “ ^ 

2 2 

t d y 

9 

r^\N^ 

N 

kn) . 

qi 


2e2/3 


pj + 



1 '- ^ 

r 13 y 

r? 


m ' 

\m) . 

Pi 


-V{ip) qj 


{N{-6 + 9e - 3e^) + M(3 - 3e)}q] + {iV(6e 
1 




^Me\{pi - qif 


4e^) - 2Me}qi{pi - qi) 
(29) 


where k is the norm of a three-dimensional wave vector since we have used the Fourier 
decomposition. We drop the subscripts I in the following and we finally sum up the two 
polarization states, or multiply results by two. Substituting background equations, we obtain 


6^C = 

4 
4 


1 

2iW 


1 e 




2e2“ 

-'-wy 






1 p^ 1 

.2ii^ ~ 2^ 
{p-qf . 


p 


(30) 


Furthermore, if we substitute the consistency relation M = (eN, this becomes 


6^C = 


jVeSa 


<^\2N^ 2e^^ J ^ V 


( ^ 2 

2CN^ 2e2"^ 


m'^M^ 


(3e 


Ce^){p-qf 


(31) 


The interaction term is proportional to {p — qY ^md do not disappear in the slow-roll limit 
where e = eo and C = 1- Thus, it is difficult to treat them analytically. Hence, we change 
the variables in the next subsection in order to have decoupled equations in the slow-roll 
limit, which allows us to calculate power-spectra analytically. 


B. Mass eigenstates 

We calculate the second order Lagrangian in the mass eigenstates by changing the vari¬ 
ables in the flavor eigenstates to those in the mass eigenstates. Then, the calculation be¬ 
comes much simpler because a massive eigenstate stemming from the interaction terms and 
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a massless eigenstate due to the general covariance do not interact with each other in de 


Sitter spacetime in bimetric gravity 


23|. 


We make a transformation of the variables {q, p) into (x, y) as 


K 





(32) 


where k = C^/^tana. We note that k is constant since ( is constant but e is time dependent 
under slow-roll approximation as we mentioned in the last part of the previous section. 
Therefore, this mixing matrix is time dependent. We can see that y corresponds to a massive 
mode since y is represented as p — q which combination appears in the mass terms and x 
corresponds to a massless mode since x is orthogonal to y. By substituting this relation into 
eq. ca. we hnd the second order Lagrangian in the mass eigenstates up to the hrst order of 
the slow-roll parameter as 


K 

4 




1-2 1-2 3rr2 

-x^ H —if H—- n 


5Cx^ + -H 


3 „2 


2e^ 


1 k‘^ / 

2 V ^ ^2 


k‘^ 


2 

ne 


-|- e 


:6Cy‘^ - 2H^^6C,xy - -ml^y"^ 


X 


1^2 


+ e 




2e 


2k^ 


-|- e- 




(33) 


in the iV = 1 gauge, where 


TUgg := m cos a- 


—(x^ + y) 


= m 


e cos^ a -I- - sin^ a 


(3 - e - (e) 


= mgg(e) -h 2s(l - eo) [3Hq - nr el cos^ a] 


(34) 


and 


tj2 2 • 2 ^0 

Hq := m sm a - . 

eo 

We have used the relation e = eH5C, in eq. fl33|) and also eq. 
we use the N = e°‘ gauge, this Lagrangian is rewritten as 


4 


52£ = 


L2 2 2 g2 


e^ o 3^,0 


(35) 
If 


2e2 


2^ ^2 + ^2- 


.2 

„2 


dC)x^ -2k 


SCx^ + T 

2 

1 


in the last line of eq. 


.2 nm ^™2 „ 2 o „,2 


6Cy^ - 2V. 


Kx'y - -ml^e^^y^ 


fi;2 + e2--- 2 


ne , 1,2 2fi;2 


;Kxy- 1 + 


fi )2 + e 2 2 V fi )2 _|_ ^2 


kw 


( 36 ) 


The primes denotes the time derivative with respect to the conformal time dehned as dt = 
e'^dp and V := a'. Furthermore, we make a scale transformation to make kinetic terms 
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canonical as 





( 37 ) 


Then, we obtain the following Lagrangian. 


= -X'2 
2 


.ly'2 ^ 1 
2 2 


2 — s + 
2 -s - 


3e2 




^eff 




1 + 


26^ 


r^C 


+ 


+2 


Ke 


/2 


+ e- 


(5C Vh"-A;" 1 + 


2k^ 








5C{-HX'Y + - rXF) 


= \x'^ 1 r2 + 3s + 

2 : 


3e^ 


r<5C 


--F'2 + 1 
2 2 L 


rj ^ 

2 ^ 3 ^ (1+1^ 




26^ 


+ e- 


r<5C 


X^ 


+ 


K^+t- 




7]^ 


-k^l + 


2k^ 


tY + e 


-M 




2Ke 




iC ( -X'Y + —XY- eXY 


(38) 

^rj rj^ / 

We have used the relations Ti' = {1—3)71 and H = "H/e" in the hrst line and 7i = —1/(1— 5)77 
in the second line. By substituting eq. fl23ll and neglecting Cl(s^), we hnally obtain 

■2 + 3dl + %r) 


fc = + i 

2 2 


--F'2 + i 
2 2 


-r 1 + 


4se^(l - eo) 


^^l + el 


2 + 3s(l 


2k§(1-£o) ^ _ (l+2g)m; 


^2 I 

KQ-l-eg 




-r 1 + 


x^ 


AskI{1 - eo)\ 


^0 + ^0 


■ 


F^ 


XF - fc^XF 


(39) 


+ X°f^-/°) (Xy + 1 

where Kq = fan a. We note that the cross terms of X and F vanish in the slow-roll limit and 
the Lagrangian becomes diagonal in the mass eigenstates. We can see that the propagation 
speed of the massless graviton differs from that of the conventional case under slow-roll 
approximation due to the difference between the lapse function of the physical metric and 
that of the other metric. The coefficients of the terms proportional to l/r]"^ also differ from 
that of the conventional case as the result of time dependence of the mixing matrix, which 
causes the modihcation of the wavenumber dependence of tensor power-spectra. 

Since we have obtained the second order Lagrangian of tensor perturbations in the mass 
eigenstates, in which the modes are decoupled from each other in the slow-roll limit, we are 
now ready to calculate the tensor power-spectra. 


11 





































IV. TENSOR POWER-SPECTRA IN BIMETRIC GRAVITY 


We calculate the tensor power-spectra generated during inflation in bimetric gravity. 
In the mass eigenstates, we can solve the free part. Then, we calculate the first order 
corrections for them by using the interaction picture and obtain tensor power-spectra in 
the mass eigenstates. Finally, we obtain the tensor power-spectra in the flavor eigenstates 
by using the relation between the variables in the flavor eigenstates and those in the mass 
eigenstates. 


A. Interaction picture 


First, we transform the Lagrangian to the Hamiltonian and quantize the system. Then, 
we calculate the correlations of the massless state and the massive state by making use of 
the interaction picture. 

The conjugate momenta of X and Y are calculated as 




dC 


Y 

X' + s- , 

V 


where we defined 


TTy 


^ = Y' 
dY' 


(40) 


4^060(1 — eo) 

s — s -- 7 ^— . 

«o + eo 

By performing Legendre transformations, we obtain the following Hamiltonian. 


(41) 


u = TlxX' + TTyW - 5‘^C 


— —X — 

2 ^ 2 


^ 2 !^^ 4 se 2 (l - eo)^ 2 + 85(1 + 


«o + eo 






1 2 1 

-TTy' -\- — 

2^2 


4sfi:^(l-eo)\ ^ V + W 


kHl + 


Kq + Cq 
2 




y2 


+ ~s(--T:xy - \XY + eXY ) . 

\ rj 7]^ / 

We note that we have defined X(k, rj) and Y (k, rj) as 

d^k 


X{x) = 


(27r)3/2 


+‘^-^X(k,r 7 ) , Y{x) = 


d^k 

(27r)3/2 


e^'^'^Y{k,r]) 


(42) 


(43) 
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We define free fields X/ and Yj where the time evolution is governed by the Hamiltonian 

4se2(l-eo)\ 2 + 3s(l + ^ 


Ho ■= [ d^k ^ 


r 1 + 


Kg + eg 


r]^ 


^7 


1 2 1 
+ 2 "’'' + 2 


Aof,2(i_.\. 2 + 3 s (1 + ^^|^) 

Pll + -^- H — 


Kg + eg 


T]- 


We separated the interaction Hamiltonian from the free one 


:= I d^ki{-Y x,Y, - ^X,Y, + k^X,Y, 


(45) 


The helds Xj and Y/ can be expanded by creation and annihilation operators 

X/(k, 7]) = ux{k, r])ak + u^ik, r])a^_^^ , (46) 


Y/(k, r]) = uyik, r])bi^ + My(/c, ni)b'_^, , 
where mode functions satisfy 


u'x + 


4se2(l-eo)\ ^ + 3s(l + ^ 


r 1 + 




Mx = 0 , 


and 


< + 




4.kS( 1 - f„)x 2 + 3s(l + 2^1^) - 


Kg + eg 


^2 


Uy = 0 


The correlation functions of these variables can be deduced as 


{Xi(k,r])Xi{k',r]')) = ux{k,r])u*x{k,r]')6^(k + k') , 


(47) 


(48) 


(49) 


(50) 


(X,(k,,,)i',(k',,,'))=0, 


(51) 


and 


{Yj{k,r])Yj{k',r]')) = uy{k,r])u*Y{k,v')d%k + k') . (52) 


The correlation functions of the original variables are calculated by using the interaction 
picture, up to the hrst order of the slow-roll parameter, as 


(X(k,7^)X(k',r;)) 

= (^/(k,??)X/(k',?7)) +i j dr7i((ihint(r7i)X/(k,r7)X/(k',r7)) - {Xi{k,rii)Xi{k!,rii)Hi^t{r]i))) 

J —00 

= {Xiik,r])Xj{k',T])) = \uxiv)f S^{k + k’) (53) 
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and 


(y(k,),)r(k',r,)> 

= {y,(k,i,)y,(k',!,)> + ; f" dr,i({ffi„,(r,i)y/(k,.))y,(k',.,))-(y,(k.r,)y,(k', 

J —OO 

= (y,(k,r,)y,(k',),)) = K(r,)|"i"(k + k') , (54) 

In the lowest order, there is no effect of interactions. However, the cross correlation gets 
corrections as 


(X(k,r/)F(k',r;)) 

= (^/(k,? 7 )h/(k',r 7 )) +i [ dr]i{{Hi,,t{r]i)Xi{k,r])Yi{k',r])) - {Xi(k,r])Yi(k',r])Hi,,t{r]i))) 


= -29 / dr]i{Hi^t{Vi)Xi{Kv)yi(^',v)) 


' —OO 

r*?? 


1 , d 


= 2s9 / dr]i / d^ki—{—Xi(ki,T]i)Yi{-ki,T]i)Xi(k,ri)Yi(k',T])) 
J-oo J Vl d7]i 


1 

vl 


--kl) (X/(ki, r]i)Yi{-ki, r 7 i)X/(k, r])Yi{k!, r])) 


= 2s5^(k + k')9 [ dr]i—Uxik,T]i)u*x{k,T])uYik,r]i)uY(k,r]) 

J-oo Vl 

— - k‘^]ux{k,r]i)u*x{k,r])uY{k,r]i)uY{k,r]) . 
^Vi ' 


(55) 


We will calculate these correlations explicitly in the long wavelength limit in the following. 


B. Power-spectra in the mass eigenstates 

We calculate the mode functions up to the order needed for calculation and obtain the 
tensor power-spectra in the mass eigenstates. We will see that the amplitude of the power- 
spectrum of the massless state remains at the end of inflation, while those of the cross 
spectrum and the power-spectrum of the massive state rapidly decay and are negligible. 


1. Power-spectrum of the massless state 

First, we calculate the power-spectrum of the massless state. We see that the power- 
spectrum looks like that of general relativity. 
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In the remote past rj —)■ —oo, the differential equations of motion for the mode functions 
become 


u'x + k\ux = 0 , u'y + kyUy = 0 , 


(56) 


where 


kx = k(^l + 


4se2(l-eo)\V2 




0 + <^0 ' 


+ ■ (57) 

V Kq tQ / 


Thus, we adopt the Bunch-Davies vacuum states as initial conditions 


ux{k,v) = 


^-ikxV 


V‘^kx 

Then, we obtain the solution of eq. fldSp as 


UY{k,r]) = 


\/2fcy 


^-ikyV 


(58) 




, - . (J<J 


sm vxx 


n=0 






2n—ux 


r(n + 1 + Vx) 


r(n + 1 - z/jc) 


(59) 


where 

3 , , 2e2(l-eo) 

4 + el 

The Hankel function (x) is dehned by 

^{1)7 5 ^yL-^.A-mz/2)^^+’^ (-1)"(z/2)2-- 

sin z/vr 

This solution reduces to 


(60) 


n=0 


r(rz, + 1 + v) r(n + 1 — z/) 


(61) 


ux{k,v) ^^—e 4 — 

2 TT 




( 62 ) 


for r] ~ 0. In the large scale limit, the power-spectrum of the massless state is given by 

d?kd?k' 


{X{z)X{z)) = 


(27r)3 
dk 47rfc^ 


e*(k+k')-z(X(k,r/)X(k',r7)) 


\ux{k,ri)f 


k (27r)^ 

dk Ank^ —7ir]T{h'x)‘^ ( —kxV\~‘^^^ 
k (2'Jiy 4 TT^ V 2 / 


( 63 ) 
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where = rj, = (z)*. Therefore, we obtain 


{x{z)x{z)) = 




:{X{z)X{z)) 


= / dilogk) 


Ho 

irMn 


-Hov) 


2s 


-kt] 


3—2ux 


V r(i) ) 


where 


A= 1 + 


4se^(l - eo)\V2 


Kn + Cn 




. -2 !^x 


(64) 


(65) 


''0 ^ '^0 

and Ho is the constant Hubble in the slow-roll limit. We note that we have used the relation 
e“ = (—obtained from the definition of s and A comes from the modification of 
the propagation speed. This amplitude is almost constant as in the conventional case, i.e. 
(xx) ~ 0{ri^) . 


2. Cross power-spectrum 


Next, we calculate the cross spectrum of the massless state and the massive state. We 
conclude the cross spectrum is negligible compared with the power-spectrum of the massless 
state. 

For calculating the cross spectrum in the first order of the slow-roll parameter, we only 
have to know the mode functions in the lowest order. The mode functions of X which we 
have calculated reduces to 


ux{k,r]) = 


\/^V kr]J 


( 66 ) 


in the slow-roll limit and the solution of eq. fH^ is written as 

y/-7rT] (1) 


UY{k,T]) = 






\/-vrr7 -2. 1 

-!-e*4 


sinh UTT ^ 

’ rn -1 


oo 

i:[« 

n=0 


,^7r 




r(n -hl + iju) 


r(n -j-1 — iju) 


sinh pvT 


E 

n=0 


^ (-l)'^(-kr]/2f^+^>^ (-lY(-kT]/2)^^-^f^ 

e 2 -e 2 


r(n -|- 1 -|- ijj) 


r(n -|- 1 — ijj) 


(67) 


where 


= 


/ 


V HI 



( 68 ) 
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Note that /i is a real number since we have proven mgg(eo) > in our previous paper |6|. 
The cross spectrum of the massless state and the massive state is 




= 2s(5''(k + k')^5 


dm 


' —OO 

rv 


= 2§6%k + k')^ dm 


1 

L^i 

-k^ 


u'xiKVi) + (\- k^)ux{KVi) 


u*x{k,r])uY{k, m)u*Y{^,^) 


■ /*oo 

■kriYi^ 


= ^(53(k + k')3f? 

dvr 


kt] 

M (—k7j\-^i^ 


e 2 


TTfl 

X < e~ 


' —kr] 

r(-ip) 


r(vp) + e 2 
f—kriyi^ _2m f—kri^^-^i^ r(vp) 

T—j (fm) 

23fJ 


ST] 


47r(m2g/hf2) 


<53(k + k') 


Svr 


]i tanh 71 ]i 


-kr]\-^i^^f3 . \ 2/. X 


+ Oiv^) 

Oiv") , (69) 


and {X{z)Y{z)) scales as t] in the leading order. Therefore, 


Similarly, {y{z)x{z)) = {{x{z)y{z)))* = 0{rf). The amplitude decays as e“^“ during infla¬ 
tion and we can neglect this amplitude compared with that of the power-spectrum of the 
massless state. 


3. Power-spectrum of the massive state 


Finally, we calculate the power-spectrum of the massive state. It turns out that it can 
also be neglected compared with the power-spectrum of the massless state. 

The mode functions of Y satisfy 


My + 


e + 


Civ) 


uy — 0 , 


(71) 


where 


Civ) := 




-2-3s 1 + 


2ft:g(l - ep) 

4 + 4 


^ — C'o + sCiiv) ■ 


(72) 


Here, Cq is dehned as Co := mlf^ieo)/H q — 2 and Ci(? 7 ) is a function of v since e has time 
dependence e = eiv) under slow-roll approximation. We define y = log(—r^), then, this 
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equation becomes 


(P d 

— -r'^Y 
d y dy 


(4e^» + C{y))uY = 0 . 


(73) 


From the definition of y, we can see y —)■ —cxo as rj ^ —0 . If we decompose the solutions 
like uy = /(?/), then we obtain 

^f+{-\ + kW" + C{y))f-0. (74) 


We can neglect the term which is proportional to because we are interested in the 
behavior of solutions around \ri\ ~ 0. We obtain WKB solutions of this equation as 


/ 


(4C-1)V4 ^"P 





since 


dCjy) 

dy 


0{s) , 


d^Cjy) 

d?y 


0{s^) . 


Finally, we obtain 


Uy 


y 

62 


{AC - 1 )V 4 


exp 





(75) 


(76) 


(77) 


up to the overall coefficient. From this, we can calculate the power-spectrum of the massive 
state 


{Y{z)Y{z)) 


d?kd^k' 


e*(k+k')'^(y(k,r7)y(k',r7)) 


(27r)3 
dkAnk^ 2 

dk Airk^ —rj 
T ( 27 r )3 {AC - 1 ) 1/2 


= 0{rj) . 


(78) 


Thus, we reach the final result 


= o(y). ( 79 ) 

As mentioned before, we can neglect this amplitude compared with that of the power- 
spectrum of the massless state as in the case of the cross spectrum. 
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C. Power-spectra in the flavor eigenstates 


We calculate the tensor spectra in the flavor eigenstates by making use of the results in 
the previous subsection. We will see all of them agree with each other and these amplitudes 
are conserved on super-horizon scales. We also see the spectral index and the amplitudes in 
the leading order. 

Using the relation between the mass eigenstates and the flavor eigenstates: 

(:;;.)(:)■ 

we can obtain the following power-spectra 

{qq) = 2 ! 2 - «e((a;|/) + (yx)) + e^yy)] , ( 81 ) 

(qp) = ^2 g2 + («Ve)(a^2/) - - K^{yy)] , (82) 

(pp) = 2 ! 2 + («Ve)((^2/) + (yx)) + (f^V^^){yy)] ■ (83) 

According to the results in the previous subsection, we And 


(qq) = (qp) = (pp) 


K 




(xx) + 0(p^ 


= / d(logk) 


K 


V ttM, 


Hn 


-Hop) 




k 


X 


f n^x) \ 
V r(|) J 


2Ho 

2 


^ «A + eS / 


6 seg(l - ep) 

4 + 4 


(84) 


The power-spectrum of the physical metric, that of the other metric and the cross spectrum 
are the same in the first order of the slow-roll parameter at the end of inflation. When we 
include both plus mode and cross mode, the results should be multiplied by two. Superfi¬ 
cially, these spectra do not seem to conserve. However, we can verify their conservation as 
follows 

dlog(g2) d r ( 2 , 2^ , f ^ 2.+ ' 

= — - logffi: + e ) + I-2s- ^ ^ 

Kq + Cq 

2 eg(l - eo) 1 


dp 


dp 


\og{K^ + e4 + (-2s^^4V^) 

k /xn 6n / 


2 ee' 


2 s- 


4 + 4 p 


H? + 

-4s'He^(l - eo) 4seg(l - Cq) 1 

p 


4 + 4 


4 + 4 


= 0 . 


(85) 
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We can read off the spectral index of the tensor power-spectrum from eg. (1841) as 


riT = —2s( 1 -|- 


2eo(l - <^o) 

nl + el 


( 86 ) 


It means that the spectra are red-tilted since eo has a value between 0 and 1 as we mentioned 
in section m The amplitudes of them in the leading order, i.e. in the slow-roll limit, are 

2 


(qq) = (qp) = (pp) = 


Kn 


4 + 4 


Ho 

TlMr, 


(87) 


The amplitude of the physical metric is suppressed compared with that in the conventional 
case since the physical metric is represented as a superposition of the massless mode and the 
massive mode and the massive mode rapidly decays. When we take the general relativity 
limit [Mf/Mg —)■ 0), the amplitude and the spectral index of the physical metric are smoothly 
connected to those in the conventional case since kq becomes inhnity. On the other hand, 
when we take the massive gravity limit {Mf/Mg —)■ cx)), the amplitude vanishes since kq 
becomes zero. 

We emphasize that the physical tensor modes and the other tensor modes are maximally 
correlated and they identically behave at the end of inflation. 


V. CONCLUSION 


The deviation of the tensor power-spectrum from the conventional case in the slow-roll 
limit is that the total amplitude is suppressed due to the mixing of the physical metric 
and the other metric. The spectral index is red-tilted compared with the conventional case 
in the leading order. The red-tilted spectrum originates in the time dependence of the 
mixing matrix, i.e. graviton oscillation. The difference between the lapse functions results 
in that the propagation speed of the massless state is larger than the speed of light and the 
amplitude is modihed in the hrst order of the slow-roll parameter. We note that the physical 
metric and the other metric are maximally correlated. This indicates that we can hnd some 
information about the other metric by observing physical quantities. 


It would be easy to extend the analysis to multi-metric gravity [2J, l25|. We have to 


calculate not only the tensor perturbations but also the scalar perturbations in order to 
connect these results to observational data such as the tensor-to-scalar ratio. Though some¬ 
one worries that this branch may suffer gradient instability during the radiation dominant 


20 








era [^, 26|, gradient instability will be avoided when we consider the bare mass is not so 
small to explain the current accelerated expansion. This situation is allowed since we have 
not only a massive graviton but also a massless graviton in the case of bimetric gravity. 
Other possibility is that we can obtain another branch which has no Higuchi ghost when 
we extend the minimal model to general models including other parameters. In those cases, 
we may need to have some extreme values for the theoretical parameters to obtain the new 
branch. We are also interested in the behavior of the tensor perturbations in the reheating 
era. It is because the other metric interacts with the scalar field in a non-trivial way through 
(. The function ( oscillates when the scaler field oscillates and this will cause the parametric 
amplihcation of the other metric. Then, the physical metric can be enhanced through the 
mixing of the physical metric and the other metric. 
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